Lifetime measurements of first excited states in ^^^^C 

H. J. Ong/'UN. Imai,2'|t| D. Suzuki/ H. Iwasaki/'El H. Sakurai/'ll T. K. Onishi/ 

M. K. Suzuki/ S. Ota/ S. Takeuchi/ T. Nakao/ Y. Togano/ Y. Kondo/'y 

N. Aoi/ H. Baba/ S. Bishop/ Y. Ichikawa/ M. Ishihara/ T. Kubo/ 

K. Kurita/ T. Motobayashi/ T. Nakamura/ T. Okumura/ and Y. Yanagisawa^ 

(^ ' ^Department of Physics, University of Tokyo, 

O '. 

C^ . Hongo 7-3-1, Bunkyo, Tokyo 113-0033, Japan 

> 

o 



c^ 



> 



o 



s 



^RIKEN Nishina Center, Hirosawa 2-1, Wako, Saitama 351-0198, Japan 



^ ■ ^ Center for Nuclear Study, University of Tokyo, 

P^ . RIKEN campus, Hirosawa 2-1, Wako, Saitama 351-0198, Japan 

"^Department of Physics, Rikkyo University, 
^ ■ Nishi-Ikebukuro 3-34-1, Toshima, Tokyo 171-8501, Japan 

^ I ^Department of Physics, Tokyo Institute of Technology, 



Ookayama 2-12-1, Meguro, Tokyo 152-8551, Japan 
(Dated: February 2, 2008) 



^ ■ Abstract 

tH- ■ The electric quadrupole transition from the first 2^ state to the ground 0^ state in ^^C was 



studied through Ufetime measurement by an upgraded recoil shadow method applied to inelastically 
scattered radioactive "'^'^C nuclei. The measured mean lifetime is 18.9 ± 0.9 (stat) it 4.4 (syst) ps, 
corresponding to a B{E2; 2f — > 0^) value of 4.3±0.2±1.0 e^fm^, or about 1.5 Weisskopf units. The 
'^ . mean lifetime of the first 2+ state in ^^C was remeasured to be 18.0 it 1.6 it 4.7 ps, about four times 

shorter than the value reported previously. The discrepancy between the two results was resolved by 
incorporating the 7-ray angular distribution measured in this work into the previous measurement. 
These transition strengths are hindered compared to the empirical transition strengths, indicating 
that the anomalous hindrance observed in ^^C persists in ^'^C. 
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I. INTRODUCTION 

Structure of neutron-rich nuclei is one of the frontiers in nuclear physics. In such nuclei, 
several exotic phenomena like a halo ll| or skin l2| structure and disappearance of the 
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traditional magic numbers [3|, U, |5|] have been discovered with the advance of experimental 
techniques and accelerators. These findings are unexpected by the standard nuclear structure 
model and precede the development of the theory. 

Recently, we reported another exotic phenomenon of extremely suppressed B{E2) value 
for the transition between the first 2^ {2f) state to the ground (Oi.) state in neutron-rich 
^^C J6|. The B{E2) was obtained by measuring the mean lifetime of the 2f state {t(2i)) 
using a new experimental technique. In general, an even-even atomic nucleus has properties 
as a quantum liquid drop. In the liquid-drop model, the B{E2) value is inversely proportional 
to the excitation energy of the 2f state {E{2f)) [7]. However, the measured B{E2) of ""^^C 
was found to deviate greatly from the value expected by the empirical formula [8i] of the 
E{2i), indicating a suppressed proton collectivity in ^^C. 

Contrary to the suppressed proton collectivity, a large neutron collectivity was suggested 
based on the measurement of the interference between the nuclear and electromagnetic 
interactions in the excitation from the ground state to the 2f state [9|. Meanwhile, the 
large quadrupole deformation length of ^^C extracted from the proton inelastic scattering 
also indicates that the neutrons contribute predominantly to the strength of the excitation 
to the 2i state, whereas the protons seem to be frozen 1G|. 

The suppressed B{E2) may indicate quenched effective charges and/or the emergence of a 
new magic number Z = 6 in the light neutron-rich carbon isotopes. For neutron-rich nuclei, 
the isovector quadrupole mode is likely to be enhanced owing to the large isospin. Such 



enhancement will give rise to quenched isovector effective charges ll|, ll2|, in particular, a 



small neutron effective charge that reduces the contribution of the valence neutron(s) to the 
B(E2) value. Indeed, the quenched effective charges have been observed in the neighboring 
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i5,i7g jiuciei |l3l.ll4l|. Moreover, in the case of the closed shell nuclei, the isoscalar polarization 



will also be reduced 



ll[ |. In this light, the emergence of the proton magic number Z = 6 



has been suggested by a shell model calculation |l5l ]. 

Besides the shell model calculation, several other microscopic models have also been 
proposed to explain the mechanism of the small B{E2) value for ^^C A calculation using 



antisyminetrized molecular dynamics (AMD) has attributed the anomalous feature to the 
opposite deformations in the proton and the neutron matters ly, [iTI]. It is interesting to 
note that calculations assuming a simple structure that ^^C is composed of ^^C +n+n have 
also reproduced the -E(2^) and the B{E2) of this nucleus 18|, ll9|, |20|. To shed light on 
the exotic phenomenon exhibited by ^^C and to scrutinize the claim for the emergence of 
the Z = 6 magic number in the neutron-rich C isotope, more experimental information 
specifically the information on the neighboring ^'^C isotope is needed. 

In this article, we report the lifetime measurement of the 2^ state in the ^^C nucleus by 
means of an upgraded recoil shadow method. The remeasurement of the lifetime of the 2^ 
state in ^^C is also reported. The r(2|') of ^^C is expected to be as long as that for ^^C since 
the energ y le vel schemes of the nuclei are almost identical. Specifically, the -^(2^^) values 
are 1766 |2l| and 1585(10) keV |22|; the one-neutron separation energies (Sn) are 4250(4) 
and 4180(30) keV, for ^^C and ^^C, respectively (211]. In addition, we have measured the 
lifetimes of the l/2~ state in ^^Be and the 3^ state in ^^N to verify the method. 



II. RECOIL SHADOW METHOD (RSM) 
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^■i'or earlier works on the recoil shadow method (RSM), we refer the readers to Ref. 23|, 
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261]. In both the previous u\ and the present work, we applied the RSM to 7 decays 



from the excited radioactive nuclei. The RSM makes use of the shadow effect of a lead shield 
placed around a reaction target on Nal(Tl) detectors of a highly segmented 7-ray detector 
array. In this method, excited nuclei are produced via inelastic scattering or fragmentation 
reaction of an intermediate-energy radioactive nuclear beam (RNB) in the reaction target. 
Each nucleus then travels a certain distance before decaying through emission of a 7 ray, 
which is detected by the Nal(Tl) detector array. The yield of the 7 rays detected by each 
Nal(Tl) detector depends on the emission point and the 7-ray angular distribution; the 
latter is governed by the alignment produced by the nuclear reaction 27|. Since the flight 
length, i.e. the emission point of the de-excitation 7 ray depends on the velocity of the 
ejectile (which is known) and the lifetime of the excited state, assuming that we know the 
angular distribution of the 7 rays, the lifetime can be determined, in principle, by observing 
the yield distribution of the 7 rays. For lifetimes of as short as a few tens of picoseconds. 



however, the yield distribution is not sensitive to the small variation of the emission point. 



Hence, determination of such lifetimes is achieved with the presence of the lead shield, which 
enhances the lifetime dependence of the yield distribution. 

In the previous work 6|], we bombarded a ^Be target with a 40-MeV/nucleon ^^C beam 
and detected the de-excitation 7 rays with only two layers of Nal(Tl) detectors placed cylin- 
drically around the beam axis. The mean lifetime was determined by comparing the ratio 
between the 7-ray yields detected in the two layers with a calculated ratio function, which 
was obtained for several mean lifetimes through simulations. However, the angular distri- 
bution of the 7 rays was not measured, and thus we had assumed the angular distribution 
obtained by a theoretical calculation. 

In this work, we have made three improvements to the RSM. First, we determined the 
lifetimes independent of the angular distribution. For the determination of each lifetime, 
two measurements of the 7-ray yields were performed: one with the lead shield installed, 
and the other without the lead shield. The 7-ray yield from the former measurement carries 
the information of both the lifetime and the 7-ray angular distribution, while that from the 
latter measurement carries mainly the information of the angular distribution. Hence, we 
can determine the lifetime by comparing the ratio of the two 7-ray yields to a simulated 
ratio function. Second, we compensated the drop (by almost two orders of magnitude) in 
the beam intensity from ^^C to ^^C by increasing the number of detectors from 32 in the 
previous setup to 130. The last improvement was achieved by using a fast ^^C beam at an 
energy of about 79 MeV/nucleon, thus increasing the flight length by a factor of about 1.4. 

As will be shown later, the results for ^^Be and ^^N demonstrate that the present RSM 
is a reliable and powerful means to determine the lifetimes of nuclear excited states with 
simple or known 7-decay schemes. Because the present RSM is independent of the 7-ray 
angular distribution, it is applicable to any type of nuclear reaction that produces excited 
nuclei. And as such, one can measure the lifetimes of excited states produced by any kind of 
reaction channels with the intermediate-energy RNB in a single experiment. In fact, we have 
also successfully measured the lifetimes of excited states in ^^C, produced by one-neutron 
knockout reaction of ^^C 1281. 



III. EXPERIMENT 

In the experiment, we measured the hfetime of the 2^ state of ^^C populated via inelastic 
scattering of the ^^C beam at 79 MeV/nucleon and those of the excited states in ^^Be and 
^^N produced through breakup reaction of the ^^C beam. We also measured the lifetime 
of the 2f state of ^^C populated via two reactions, namely the inelastic scattering of a 
72-MeV/nucleon ^^C beam and the breakup reaction of the 79-MeV/nucleon ^'^C beam. 
Furthermore, we measured the angular distribution of the 7 rays for the 2^ state of ^^C 
inelastically excited at 40 MeV/nucleon in order to incorporate it into the experiment data 
in Ref. 6|. 

The experiment was performed at the RIKEN accelerator research facility. Secondary 
beams of ^^'^^C were produced in two separate measurements through projectile fragmenta- 



tion of an 110-MeV/nucleon ^^Ne primary beam, and separated by the RIPS beam line 29|. 
Particle identification of the secondary beam was performed event-by-event by means of 
the time-of-fiight (T0F)-Ai5 method using two 1.0-mm-thick plastic scintillation counters 
located at the second and final focal planes of RIPS. The beams were directed at a 370- 
mg/cm^ ^Be target placed at the exit of the RIPS beam line. Two sets of parallel plate 
avalanche counters (PPACs) were placed upstream of the target to measure the position 
and angle of the projectile incident upon the target. The ^^.isq beam had typical intensities 
of 6.5 X 10"^ and 2.3 x 10^ particles per second, respectively. 

Outgoing particles from the target were identified by the AE-E-TOF method using a 
plastic scintillator hodoscope |30J located 3.8 m downstream of the target. The hodoscope, 
with an active area of 1.0x1.0 m^, consisted of thirteen vertical Aii^-scintillator slats and 
sixteen horizontal ii^-scintillator bars with 5.0 mm and 60.0 mm thicknesses, respectively. 
The scattering angles were determined by combining the hit position information on the 
hodoscope with the incident angles and hit positions on the target obtained by the PPACs. 

In order to implement the RSM concept, a thick 7-ray shield was placed around the 
target, as shown in Fig. [T] The shield was a 5 cm-thick lead block with an outer frame of 
24x24 cm^ and an inner hole of 5.4 cm in diameter. The inner hole surrounded the beam 
tube housing the ^Be target. We removed the lead shield during measurement of the angular 
distribution of 7 rays. For the sake of later discussion, the 2;-axis is defined as the beam 
direction, which is close to the flight direction of the de-exciting nucleus. The origin of the 
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FIG. 1: Schematic view of 7-ray detectors. A beryllium target is surrounded by a 5 cm-thick lead 
shield, and 10 layers of Nal(Tl) scintillators (LI ~ LIO) placed cylindrically around the beam axis. 
For clarity, only part of the detectors and lead shield is shown in the inset. 

2;-axis, z = 0.00 cm, was taken at the upstream edge of the lead shield. We have carried out 
both measurements, namely the measurements with and without the lead shield, with the 
target placed at 2 = —0.15 cm and z = 2.05 cm. We refer to these setups as the "upstream 
setup" and the "center setup" hereinafter. 

The 7-rays from the excited nuclei in-flight were detected by an array of 130 Nal(Tl) 
detectors, which form part of the DALIl 31| and the DALI2 [32]. The Nal(Tl) detectors, all 
of which are rectangular in shape, are of three different sizes: 4.5x8x16 cm^, 4x8x16 cm^, 
6.1 X 6.1 X 12.2 cm'^. The array was divided into 10 layers, labeled Ll-10 as shown in Fig. [U 
with each layer consisting of 10 ~ 18 detectors arranged coaxially with respect to the beam 
direction. The layers are closely packed to cover polar angular ranges of 15° ~ 85° and 
100° ~ 170° in the laboratory frame. The detectors were mounted on and supported by nine 



3-nim-thick aluminuin plates. The distances of the center of the detectors from the center 
of the Be target ranged from 28 cm to 56 cm. 

In the present work, we counted the number of full-energy-peak events detected by each 
layer during the measurements with and without the lead shield. For convenience, the 
numbers obtained with the i-th layer during the former and the latter measurements are 
denoted by iV^Pb ^^^ ^woPb- The iV^Pb ^^^ ^woPb were obtained by fitting the measured 
7-ray energy spectra with response functions obtained with simulations, plus 7-ray spectra 
for the ^^C isotope as backgrounds. The spectra for ^^C were selected because (i) all excited 
states in ^^C lie at energies above 6 MeV, and there is no significant 7 line around 1.6 MeV, 
(ii) the ^"^C isotope was produced through projectile fragmentation reaction of the ^^C (or 
^^C) in the secondary target, and as such its "background spectrum" resembles that of for 
i^C (or i^C). 

As mentioned earlier, a key point of the present work is the elimination of the dependence 
of the measured lifetime on the 7-ray angular distribution. For this purpose, we determined 
the deficiency of the 7-ray yields due to the lead shield. The deficiency (D*) of the i-th layer 
detectors is defined as the ratio between the yields detected with and without the shield, i.e. 

D' = /.iv;pb/A^;opb, (1) 

where /;, is the normalization factor for different total number of beam particles in the 
two measurements. The lifetime was determined by comparison of the measured deficiency 
(-^exp) with the simulated one (Dl:^^). 

IV. MONTE CARLO SIMULATIONS 

The simulated deficiency of each layer as a function of various lifetime for the respective 
de-excitation 7 rays was obtained by performing Monte Carlo simulations using the GEANT 
code 



33|. We have taken into account the geometry of the experimental setup, including 
the shape of the detectors, in the simulations. The geometry was checked by performing 
separate measurements in which ^^"^Cs, ^^Na and ^°Co standard sources emitting 662-keV, 
1275-keV, 1173-keV, and 1333-keV 7 rays were placed at several positions from z = —0.15 
to 5.15 cm. The deficiencies of all layers measured for the respective target position were 
reproduced by the simulation within accuracies of ± 3% and ± 7% for layers with -D* — 0-2 
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and Dgxp < 0.2, respectively. 

The simulation was then applied to the case of 7-rays emitted from the de-exciting parti- 
cles in flight. For this simulation, we have considered the experimentally obtained parameters 
such as the energy and emittance of the projectile, the angular spread due to inelastic scat- 
tering and multiple scattering, and the energy loss in the target. Apart from determining 
the -Dgjjjj, the response functions obtained from the simulation were also used for the fitting 
of the experimental 7-ray energy spectra. Note that the strong dependence of Dl^^ on the 
lifetime is exhibited in Fig. [3l Fig. [5l Fig. [7] and Fig. [91 

V. RESULTS 

We deduced the -Dgxp using Eq. [1] and determined the lifetimes for the respective excited 
states. In the case of the inelastic scattering with the 40-MeV/nucleon ^^C beam, we revised 
the lifetime of the 2^ state reported previously [6] by incorporating the measured 7-ray 
angular distribution. The results from these measurements are presented. 

A. Lifetime of 2| state in ^^C 

We show in Fig. [2] the 7-ray energy spectrum measured with the L6 Nal(Tl) detectors, 
which were located at polar angle of about 90° in the center of mass frame. The spectrum 
was obtained with the center setup and without the lead shield. The significant peak around 
1500 keV corresponds to the 1585-keV line from the 2]*" -^ Ogg transition, while the minor 
peak around 1000 keV corresponds to the 919-keV line from the transition between the 
2504-keV excited state and the 2f state. The level scheme for ^'^C as shown in the inset 



of Fig. [2] has been proposed recently based on the in-beam 7-spectroscopy 22]. The two 7 
lines observed correspond to the transitions shown by the bold arrows in the level scheme. 
The two significant peaks around 200 keV and 300 keV are from the known transitions in 
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34| : the 7-ray energies shown are taken from Ref. 28 1. 

By considering the sum spectrum of all detectors, the cascade contribution from the 

second excited state to the 2f state was determined to be about 9% relative to the 2^ — > Ogg 

transition. Since the error that might arise from the cascade contribution is negligibly 

small compared to the systematic error mentioned below, we have neglected the cascade 
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FIG. 2: 7-ray energy spectrum for the inelastic scattering of the 79-MeV/nucleon C nuclei on 



^Be. The 7 rays were detected by the L6 Nal 
scheme of ^^C and the known 7-ray transitions 
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detectors. The inset shows the energy level 



contribution in our analysis. 

The simulated deficiencies were obtained for several assumed mean lifetimes ranging from 
r(2j^) = ps to 40 ps. As an example, we compare the D^^^ value with the -0*;^ values 
obtained for r(2+) = 0, 20, and 40 ps. The D^^p's (filled circles) and the D^^'s for r(2+) = 
ps (dashed line), 20 ps (solid line) and 40 ps (dotted line) are plotted in Fig. [3] for the (a) 
center and (b) upstream setups. Note that we have omitted the experiment data for the 



layers with A^^ 



wPb 



0, which do not contribute to the determination of the mean lifetime. 



The figures indicate that the mean lifetime of the 2^ state locates around 20 ps. The mean 
lifetime was determined by searching the x^ minimum; the x^i"^) is given by 



x\r) = y: 
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FIG. 3: The Dl^^s, denoted by the filled circles, of the respective layers for the 2f — *• Oi, transition 
in ^^C as compared with the simulated values. The dashed, solid and dotted lines represent the 
^sim values calculated for T{2f) = 0, 20, and 40 ps, respectively, (a) The D^ distribution for the 
center setup, (b) The D* distribution for the upstream setup. The insets of the figures show the 
reduced x^ distributions as functions of r(2^). 

where the SDl represents the statistical error. The x^ distributions for t{2^) from ps 
to 40 ps, as shown in the insets of Fig. [3l give t{2^) = 19.9 ±1.0 ps and 14.9 ±2.1 ps for 
the center and upstream setups, respectively. We adopted the weighted mean of these two 
values, 18.9 ± 0.9 ps, as the mean lifetime for the 2'^ state in ^^C 

The systematic error mainly composed of two uncertainties. One is attributed to the 
discrepancy between the measured deficiency and the simulated one, and the other is due to 
the uncertainty of the target position, which was about 0.5 mm. The former was estimated 
to be 0.5 ps by changing the DI^j^{t) randomly within ±3% (or ±7%) and observing the 
deviation of the lifetime thus obtained. For the latter, the uncertainty in the target position, 
which corresponds to the uncertainty in the emission point, resulted in an uncertainty of 
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about 4.4 ps for the ^^C ejectiles traveling at about 38% of the speed of hght. The resultant 
systematic error, defined as the root sum square of these two values, was determined to be 
4.4 ps. The mean lifetime thus obtained is 18.9 ± 0.9(stat) ± 4.4(syst) ps. For simplicity, 
the notations for the statistical (stat) and the systematic (syst) errors will be omitted from 
now onwards. 

B. Lifetime of 1/2 state in ^^Be 

Figure m shows an example of the 7-ray energy spectrum measured with the L6 Nal(Tl) 
detectors in coincidence with the ^^Be ejectiles. The spectrum was obtained with the center 
setup and without the lead shield. Only one peak that corresponds to the transition of 
the 1/2^ -^ 1/2+ is observed at around 300 keV. The 1/2^ state is the only known bound 
excited state that decays through 7 transition as shown by the level scheme [21| in the inset 



of Fig. m The mean lifetime is r(l/2^) =0.166(14) ps 2l|], which is much shorter than the 
mean lifetime of ^^C{2'f). By applying the RSM to this short lifetime, we verified the lower 
limit of the dynamic range of the method. 

The -Dexp ^^d the Dl:^^ values for r(l/2^) = 0, 15, ad 30 ps are shown in Fig. [5] for 
the two different target positions. As is obvious from the plots, the deviation of the Dl^^ 
from the -D* becomes larger with longer assumed mean lifetime. The simulated deficiencies 
for r(l/2~) ranging from ps to 30 ps were used to obtain the x^ value as a function of 
r(l/2^). From the x^ distribution, the mean lifetimes were determined to be 3.5 ± 1.3 ps 
and 9.5 ± 2.8 ps for the center and the upstream setups, respectively. The weighted mean 
of the two values, i.e. 4.6 ±1.1 ps was adopted. 

The sources of the systematic error are similar to those of the case of ^^C Nonetheless, 
it is sufficient to consider only the one due to the uncertainty in the target position of 
±0.5 mm, which was dominant. For the ^^Be ejectiles with velocity about 37% of the speed 
of light, the associated uncertainty in the determined mean lifetime is about 4.5 ps. Thus, 
the resultant mean lifetime for the l/2~ state is 4.6 ± 1.1 ±4.5 ps. This result is consistent 



with the reference value of 0.166(14) ps 2l|. The large error shows that the RSM is not 



suitable for determination of the lifetime below 10 ps. 
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FIG. 4: 7-ray energy spectrum in coincidence with the ^^Be ejectiles. The 7 rays were measured 
with the L6 Nal(Tl) detectors. The inset in the figure presents the energy level scheme and the 
lifetime of the first excited state of ^^Be. 

C. Lifetime of 3" state in ^^N 



The validity of the RSM was also tested by measuring the known lifetime of the 3 
state of ^^N. The excited 3~ state decays to the ground 2~ state with the mean lifetime 



2l[ |. The Doppler-corrected 7-ray energy spectrum measured by all 



of r(3-) = 132(2) ps 

Nal(Tl) detectors is shown in Fig. [61 A minor and a prominent peaks are observed around 
400 keV and 300 keV, respectively. The 400-keV peak corresponds to the 397-keV 7 line 
from the 1~ -^ 2^ transition, while the 300-keV peak includes two 7 lines of 277 keV and 
298 keV from the 1~ -^ 0^ and 3^ -^ 2^ transitions, as shown by the level scheme 21| in 
the inset of Fig. [61 

As is evident from Fig. [6l the 277-keV and the 298-keV peaks were not resolved in 
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FIG. 5: The Z^g^p's, denoted by the filled circles, of the respective layers for the l/2~ -^ 1/2+ 
transition in ^^Be as compared with the simulated values. The solid, dashed and dotted lines 
represent the Dl.^^ calculated for r(l/2~) = 0, 15, and 30 ps, respectively, (a) The D^ distribution 
for the center setup, (b) The D* distribution for the upstream setup. The insets show the reduced 
X^ distributions as functions of t(1/2~). 



the energy spectrum. Hence, to determine the iV^opb ^^"^ ^4pb ^^^ ^^^ 298-keV peak, it 
was necessary to determine the contribution of the 277-keV line in the 300-keV peak. We 
determined the contribution by fitting the energy spectrum of all Nal(Tl) detectors with the 
simulated response function for the 298-keV 7 line, and a combined response function for 
the 277-keV and 397-keV lines, taking into account the branching ratio 2l|] of the decays 
from the 1~ state. The contribution of the 277-keV line was determined to be about 25%. 
By fixing the 277-keV contribution to 25%, we fitted the energy spectrum for each layer with 
the simulated response functions and the 7-ray spectrum for the ^^C isotope as background. 
Figure [7] shows the Dl^^ and the Dl:^^ values simulated for r(3~) = 100, 140 and 180 ps for 
the two target positions. The plots clearly indicate that the lifetime locates around 140 ps. 
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FIG. 6: Doppler corrected 7-ray energy spectrum in coincidence with the ^^N ejectiles. The 7 
rays were measured with ah Nal(Tl) detectors. The inset shows the energy level scheme and the 
known mean lifetimes of the exited states in ^^N. 



The simulated deficiencies for t(3~) ranging from 115 ps to 175 ps were used to obtain tlie 
X^ values as a function of t(3^). From the x^ distributions, we obtained t(3^) = 137 ±4 ps 
and 136 ± 3 ps for the center and the upstream setups, respectively. The weighted mean of 
the two values thus obtained is 136 ± 3 ps. 

Unlike the cases for ^^C and ^^Be, the systematic error is mainly due to the uncertainty 
in the determination of the 277-keV contribution in the 300-keV peak. This uncertainty, 
determined to be about 5%, is ascribed to the uncertainty in the energy calibration of 
the Nal(Tl) detectors, which was about 2 keV at 300 keV. The resultant systematic error 
distributes from —7 ps to +11 ps. Hence, the mean lifetime for the 3~ state was determined 
to be 136 ± 3^7^ ps, which is in good agreement with the reference value of 132(2) ps 21|. 

We shall note that the energies of the 7 rays emitted by ^^Be and ^^N are well below 
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FIG. 7: The D* 's, denoted by the filled circles, of the respective layers for the 3~ — > 2~ transition 
in ^^N as compared with the simulated values. The dashed, solid and dotted lines represent the 
^sim calculated for t(3^) = 100, 140, and 180 ps, respectively, (a) The D^ distribution for the 
center setup, (b) The D^ distribution for the upstream setup. The insets show the reduced x^ 
distributions as functions of t(3~). 



that of from the 2^ state in ^^C, which is in the region of 1300 - 2300 keV depending on 
the layer. To rule out any possible energy dependence, we have also measured the mean 
lifetime of the 2^ state in ^^B. The result is consistent with the value of 260 fs 2l(] given in 
the literature. 



16 



72 MeV/nucIeon 







1000 2000 3000 4000 1000 2000 3000 4000 

E],'^ [keV] 



FIG. 8: 7-ray energy spectra for ^^C obtained with the L6 Nal(Tl) detectors during the mea- 
surement with the center setup and without the lead shield, (a) The inelastic scattering of ^^C at 
72 MeV/nucleon. (b) The breakup reaction of ^'^C to ^^C at 79 MeV/nucleon. (c) The inelastic 
scattering of ^^C at 40 MeV/nucleon. The level scheme with the known transitions in ^^C is shown 
in the top right panel. 

D. Lifetime of 2^ state in ^^C revisited 



1. Inelastic scattering of 72- MeV/nucleon C beam 



Figure [8] (a) shows an example of the 7-ray energy spectrum in coincidence with the ^^C 
ejectiles obtained with the L6 Nal(Tl) detectors during the measurement with the center 
setup and without the lead shield. The significant peak around 1800 keV corresponds to 
the 1766-keV 7 line from the 2^ — > Oi, transition. The two peaks around 300 keV and 



17 



the peak around 740 keV are from the known transitions in ^^C 22|, |28|, |3J] and ^^C 21| . 
respectively. The minor peak around 1000 keV, which was also observed at the same energy 
in the spectra of the other layers as well as in the spectra of other reaction channels, is likely 
to correspond to the 980.7-keV [2l| line from ^Li produced through fragmentation of the ^Be 
reaction target. Although no notable peak can be observed around 2300 keV in the figure, 
the sum spectrum of all Nal(Tl) detectors exhibits a small peak, which corresponds to the 



cascade transition from the higher excited state(s) as shown by the energy level scheme 21 1 
in the top right panel of Fig. [HI These higher excited states cannot be identified due to 
the resolution of the Nal(Tl) detectors. Nonetheless, by fitting the sum spectrum with a 
simulated response function for 2300-keV 7 rays, the cascade contribution was determined 
to be 9% relative to the 2'f -^ Ogg transition. Similar to the case for ^^C, this cascade 
contribution is negligibly small, and thus was not considered in our analysis. 

The experimental deficiency of each layer is compared with the simulated ones of 
r(2j^) = 0, 20, and 40 ps in Fig. [9] for the two target positions. The simulated curves 
ranging from t(2^) = ps to 40 ps were used to obtain the x^ values as a function of 
t{2i). The x^ distributions give t(2^) = 18.3 ± 1.8 ps and 14.8 ± 3.9 ps for the center and 
the upstream setups, respectively. Taking the weighted mean of the two values, the mean 
lifetime was determined to be 17.7 ± 1.6 ps. This value is about four times shorter than the 
value reported previously [6|. 

The sources of the systematic error are similar to those for ^^C Based on the same 
analysis, the systematic errors attributed to the uncertainties of the -Dg^j^ and the uncertainty 
of the target position were determined to be 0.5 ps and 4.6 ps, respectively. Adopting the 
root sum square of these two values as the systematic error, the mean lifetime thus obtained 
is 17.6 ± 1.6 ± 4.6 ps. 

2. Breakup reaction of 79-MeV/nucleon ^^C beam 

Figure [8] (b) shows the 7 ray energy spectrum measured by the L6 Nal(Tl) detectors. 
The spectrum was obtained with the center setup and without the lead shield. A significant 
peak around 1800 keV corresponds to the 2^ — > Oi, transition of 1766 keV. Compared 
with the energy spectrum of the inelastic scattering, the cascade transition observed around 
2300 keV was enhanced. By considering the sum spectrum for all Nal(Tl) detectors, the 
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FIG. 9: Results for the inelastic scattering of 72-MeV/nucleon ^^C beam. The Z)* 's, denoted by 
the filled circles, of the respective layers for the 2^ -^ 0^^ transition in ^^C as compared with the 
simulated values. The dashed, solid and dotted lines represent the Dl^^ calculated for t{2'^) = 0, 
20, and 40 ps, respectively, (a) The D^ distribution for the center setup, (b) The D* distribution 
for the upstream setup. The insets show the reduced x^ distributions as functions of t{2'^). 

cascade contribution was determined to be 22% relative to the 2^ -^ Oi, transition. 

In contrast to the case for inelastic channel, the cascade contribution was quite large. 
A possible sizable lifetime for the higher excited state (s) might affect the outcome of the 
lifetime determination for the 2^ state. Hence, to determine the t(2^) and simultaneously 
take into account the cascade contribution, we carried out a minimization two-parameter 
X^ analysis. The x^ values were obtained with the measured -Dgxp's and the D^jj^'s. Here, 
the -Dgim's were obtained for T{2f) and the lifetime of the higher excited state, denoted by 
r(cascade) hereinafter, from ps to 40 ps. A minimum was observed around t(2^) = 20 
ps and r(cascade) = ps. The mean lifetimes for the 2f state were determined to be 
20.2 ± 8.7 ps and 16.9 ± 16.8 ps for the center and the upstream setups, respectively. Taking 
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the weighted mean of the two values, we obtained t(2^) = 19.5 ± 7.7 ps. 

The sources of the systematic error are similar to the ones for the inelastic channel. The 
error attributed to the uncertainties of the -Dgim's is 0.5 ps. Meanwhile, the speed of the 
^^C nuclei was about 37% of the speed of light. Hence, the systematic error ascribed to the 
uncertainty in the target position was about 4.5 ps. Taking these two systematic errors into 
consideration, the resultant mean lifetime for the breakup channel becomes 19.5±7.7±4.5 ps. 
This value is consistent with the value determined from the inelastic channel. 

3. Inelastic scattering of 40-MeV/nucleon ^^C beam 

The de-excitation 7 rays from the ^^C nuclei inelastically excited with the Be target were 
measured at 40 MeV/nucleon. Figure [8] (c) shows the 7-ray energy spectrum measured with 
the L6 Nal(Tl) detectors. The spectrum was obtained with the upstream setup and without 
the lead shield. Since no notable peak was observed around 2300 keV, we have neglected 
the effect of the cascade contribution on the angular distribution of the 2t -^ Oi, transition. 

In the previous experiment [6|], we determined the t{2^) by comparing the relative yields 
of de-excitation 7 rays measured with two layers of Nal(Tl) detectors with a simulation 
assuming an angular distribution of 7 rays. The layers, labeled Rl and R2, were located 
just upstream of the target with central angles of 135° and 116° in the center of mass frame, 
respectively. 

The measured angular distribution of 7 rays emitted from the 2^ state is shown in 
Fig. [TOl together with the calculated distributions obtained with the ECIS79 [35[ code using 



two optical-model parameter sets, OMl [36|] and 0M2 [37|. These calculated distributions 
were used to determine the t{2^) in the previous work. From the figure, it is obvious that 
the calculations fail to reproduce the experimental data, especially at around 90°. In order 
to determine the angular distribution quantitatively, we fitted the data with the following 
function: 

W{e) = l/(47r)(l + aP2{cose) + bP^icose)). (3) 

The -Pi;(i=2,4) in the above function are the Legendre polynomials, while the a and b are the 
coefficients. The best-fitted result is shown in Fig. [10] by the solid curve. 

The measured 7-ray distribution gives iy(135°)/ iy(116°) = 1.00 ± 0.01. On the other 
hand, the calculated distributions yield W {135°) /W {116°) = 0.91 for OMl and 0.88 for 
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FIG. 10: The angular distribution of the de-excitation 7 rays from the 2^ state in ^^C. The 
sohd hne indicates the best fitted angular distribution. The dashed and dotted lines represent the 
calculated distributions for OMl and 0M2, respectively. 

0M2, which made the simulated R1/R2 ratio smaller. As a result, a longer mean lifetime 
was deduced. Figure [TT] shows the simulated R1/R2 ratio as a function of the t{2i) by 
incorporating the experimental angular distribution. The dashed lines represent the original 
simulated R1/R2 ratio. The R1/R2 ratio measured in the previous work for the two target 
positions, z = 0.0 and 1.0 cm, are shown by the hatched zones. The overlapped region 
between the experimental R1/R2 ratio and the simulated lines corresponds to 51 ± 21 ps 
and 20 ± 18 ps for z = 0.0, 1.0 cm, respectively. The resultant mean lifetime of 34 ± 14 
was obtained by taking the weighted mean of these two values. The systematic uncertainty, 
attributed mainly to the geometrical uncertainty, was estimated to be 20%. Hence, the 
revised lifetime is 34 ± 14 ± 7 ps. This value is shorter than the previously reported value of 
77 ± 14 ± 19 ps [6], but is consistent with the other two experimental results given above. 



VI. DISCUSSION 



The obtained t{2^) for i6,i8q g^j.^ summarized in Table [H In the case of ^^C, all three 
t(2'1) values obtained from the three measurements are in consistent with each other. We 
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FIG. 11: Two solid lines represent t(2^) vs. R1/R2 curves obtained by Monte Carlo simulation 
including the experimental angular distribution of de-excitation 7 ray for target positions of z = 0.0 
and 1.0 cm. The dashed lines represent the original simulated R1/R2 ratio [6]. The hatched zones 
represent the experimentally determined R1/R2 ratios ^] for the two target positions. 

TABLE I: Summary of the mean lifetimes of the 2f states in ^^C and ^^C, and the corresponding 
B{E2) values. The superscripts a, 5 and c indicate inelastic channel at 72 MeV/nucleon, breakup 
channel at 79 MeV/nucleon, and inelastic channel at 40 MeV/nucleon. 



r(2+) [ps] 



B{E2) [eHm 



2f^4l 



18c 
16(~;c 



18.9 ±0.9 ±4.4 4.3 ±0.2 ±1.0 

17.7 ±1.6 ±4.6 2.7 ±0.2 ±0.7 

19.5 ±7.7 ±4.5 2.4 ±0.9 ±0.6 

34 ±14 ±7 1.4 ±0.6 ±0.3 



adopt the weighted mean of the three values by considering only the statistical errors. For 
the systematic errors, since they are comparable in all three cases, we adopt the largest 
value (26%) as the resultant systematic error. The t(2^) value for ^^C thus obtained is 
18.0 ± 1.6 ± 4.7 ps. 

The T{2f)'s determined in the present work correspond to B{E2) values of 2.6 ± 0.2 ± 
0.7 e^fm^ and 4.3 ± 0.2 ± 1.0 e^fm^ for ^^'^^C, respectively. Although the B{E2) value for 
^'^C is almost twice as large as that of ^^C, it is comparable to the B{E2) = 3.7 e^fm'^ of 
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FIG. 12: (a) B{E2) values in W.u. and (b) ratios between the experimental B{E2) values and 
B{E2)sys calculated using the empirical formula |g] for even-even nuclei with A < 50. The filled 
circles denote the values of ^^'^^C, and the open circles represent data for other open-shell nuclei. 
The open squares and open diamonds denote the proton- and neutron-closed-shell nuclei, while the 
filled diamonds represent the double magic nuclei. The dashed lines are intended as eye guides. 



the closed-shell ^^C nucleus. Both the energy of the 2^ state and the B{E2) value remain 
small in ^^C, clearly indicating that the phenomenon of hindered E2 strength observed in 
""^^C persists in ^^C 

In Fig. [12] (a), the B{E2) values obtained for i6,i8q ^^^^ compared with all known B{E2) 



values for the even-even nuclei with A < 50 



38|. Nuclei with open shells tend to have B{E2) 



values greater than 10 W.u., whereas nuclei with neutron- or /and proton shell closure tend to 
have distinctly smaller B{E2) values. Typical examples of the latter category are the doubly 
magic nuclei, ^^O and ^^Ca, for which the B{E2) values are known to be 3.17 and 1.58 W.u., 
respectively. The present B{E2) values for i^.isq ^^^^ -^-^ ^^^^ ]^ g W.u., respectively, which 
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are even more suppressed than those of the doubly magic nuclei although they are supposed 
to be open-shell nuclei. 

The strong hindrance of the i^'i'^C transition can also be illustrated through comparison 

11 
with an empirical formula based on a liquid-drop model [8]. The empirical formula can be 

expressed by 

S(E2)sys = 6.47 X Z^A-^-^^E{2+)-\ 

The experimental B{E2) values relative to B{E2)sys are plotted in Fig. [TWb). As noted in 



Ref. [38J, the B{E2)/B{E2)sys ratios for most of the open-shell nuclei fall around 1.0, being 



confined between 0.5 and 2.0. Even for the closed-shell nuclei, the ratio remain larger than 
0.20. Thus, the ratios of 0.14 and 0.21 for i6,i8q ^^^^ exceptionally small. In particular, the 
ratio for ^^C with E{2'^) = 7012 keV 2l|] is as large as 0.68, suggesting different mechanisms 
for the small B{E2) in ^^C and the suppression of the B{E2) values in ^^'^^C As in the case 
of ^^C [qI, [10(], the observation of the small B{E2) value in ^^C despite the lowering of the 
E(2^) may imply a neutron-dominant quadrupole collectivity in ^^C 

Regarding the suppressed B{E2) values for ^^'^^C, calculations have been performed in 
several theoretical frameworks. Figure [13] shows the experimental and some of the theoretical 
B(E2) values for carbon isotopes from ^''C to ^°C. 

n n 

Calculations using the antisymmetrized molecular dynamics (AMD) [161 . Il7l | and the 



multi-Slater determinant AMD (AMD+MSD) 39j have reproduced the trend of the small 
B{E2) values in ^^'^^C In the framework of the AMD calculation, the presence of opposite 
deformations in the neutron and proton are said to be accounted for the small B{E2) in 
the i^'^'^C isotopes. Meanwhile, another calculation using the deformed Hartree-Fock wave 



function 40|] has reproduced the B{E2) values in i6,i8q rather well, although it fails to 



reproduce the one in the neutron-closed-shell nucleus ^^C. 



The shell model calculations [15|, |4l| have also predicted suppressed B{E2) values for 
^^'^^C In the case of the calculation using a Hamiltonian with the strengthened tensor 
interaction jl5|, the energy gap between (vr0pi/2)-(vr0p3/2) becomes as large as the gap for 
(7rlsi/2)-(vr0pi/2), when the neutron number is 8 and 10. This large gap reduces the proton 
matrix element. As a result, the Oi. — *> 2^ transition is mainly contributed by the neutrons. 
On the other hand, the dominant configuration of the valence neutrons in the sd-shell is lsi/2, 
and since the valence neutrons in the s-orbit spread out widely, the effective charges become 



small. In fact, the effective charge of neutron has been calculated to be e„ = 0.2e [13], which 
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is smaller than the standard value e„ = 0.5e used for the sd-shell nuclei. The other shell 
model calculation of the no-core type has reproduced successfully the present B{E2) values 
for i6,i8q ^j^]^ g^ small neutron effective charge of 0.164e 4l|. In the calculation, a quenched 
proton transition is expected. The occupation numbers of proton in the Op3/2 orbits change 
by less than 1% for the ground and the 2^ states in ^^C |4l!i] and ^^C 42|. Thus, these two 
calculations give the same picture for the suppressed B[E2) value. Note that the former 
shell model calculation predicts that the energy gap of {ttpi/2)-{'^P3/2) decreases from ^^C 
to ^^C As a result, proton matrix element will be enhanced, leading to larger B[E2) values 
for 20.22c. 

While the above theoretical models offer possible interpretation for the hindered E2 
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transitions observed in ^^'^^C, more theoretical and experimental studies on the ground 
states and other excited states, e.g. 2^ state in ^°C, are necessary for a more unified and 
complete understanding of the structure of the neutron-rich carbon isotopes. 

Finally, we comment on the experimental results of ^^C reported in Ref. |9|, [lO|. A smaller 



■•B- 



The underestimation 



B{E2) value consistent with the one in Ref. ^] was reported in Ref 
of the B(E2) value in Ref. 9(] may indicate the need for a microscopic approach in analy zing 
the reaction data. In fact, a subsequent analysis 43|] using the AMD wave functions 17 1 
in the microscopic coupled-channels calculations has indicated a larger B{E2) value of 1.9 
e^fm^. As for the work on the inelastic proton scattering on ^^C 10|, the quadrupole 



deformation length was extracted from the experiment data. Using this deformation length 
and the B{E2) value in Ref. |6|, the ratio of the neutron and proton quadrupole matrix 
elements {Mn/Mp)/{N/Z) was determined to be 4.0 ± 0.8. With the present B{E2) value, 
the {Mn / Mp) / {N / Z) value for ^^C becomes 1.9±0.4, which is still very large and comparable 
to the value of ^'^O [44]. In fact, the M„ value of about 11 fm^ deduced for ^^C remains 
much larger than the Mp (= ^JbB{E2)/e^) value, which was deduced to be about 4 fm^ 
using the present B{E2) value. 

VII. SUMMARY 

The lifetime of the 2f state in ^^.isq -were successfully measured using the upgraded RSM 
with 10-layer Nal(Tl) array, incorporating the inelastic scattering and breakup reaction at 
around 75 MeV/nucleon. The 7-ray angular distribution for the inelastic scattering of ^^C 
at 40 MeV/nucleon was also measured. Incorporating this angular distribution into the 
measurement reported previously [Q], the t{2i) of ^^C was found to be in consistent with 
the present result. The T{2f) values for ^^.isq \^\i^^ determined were as long as around 20 ps, 
indicating that the anomalous suppression of B{E2) observed in ^^C persists in ^^C In the 
framework of shell model calculation, the suppressed B{E2) values can be attributed to 
the small effective charges and the widening of the energy gap between the vr(pi/2)-vr(p3/2) 
orbitals. The present results, together with the small B{E2) values for ^^C, suggests a 
possible proton-shell closure in the neutron-rich 14,16,18q nuclei. 
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